Effect of inplane electric field on magnetotransport in helical metal 
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The existence of helical surface states in a bulk insulator, with anomalous magneto-electric prop- 
erties, is a remarkable new development in solid state physics. However clear signatures in transport 
measurements are lacking. In this paper we report on novel phenomena which provide a route to de- 
tecting massless Dirac fermions. The linear dispersion of the the fermions leads to a form of Lorentz 
invariance, with the Fermi velocity playing the role of the velocity of light (c). In a crossed elec- 
tromagnetic field, the single particle states form Landau levels, whose wave-functions and energies 
depend on the applied in-plane electric field but not their degeneracies. We predict the existence 
of oscillations of conductivity as a function of both the chemical potentials and applied electric 
field. This is an unique signature of the linear dispersing states due to an even-odd effect where the 
conductivities of the odd Landau levels are larger than the neighboring even ones. The electric field 
dependence also leads to testable signatures in steady state thermopower of the non-trivial states. 



Helical Dirac fermions, massless relativistic charged 
particles with spin locked to their linear momentum, ex- 
ist on the surface of three dimensional topological insu- 
lators (TI)0, 0] and experimentally observed^ by an- 
gle resolved photoemission spectroscopy. Novel phenom- 
ena such as the existence of Majorana fermions in the 
presence of superconductor-ferromagnet interface on the 
surface [J], and a realization of a magnetic switch by 
tuning the conductivity with a proximate ferromagnetic 
film [a], have been conjectured to be supported in these 
systems. The strong coupling of spin and orbital degrees 
of freedom is the source of these proposals. Magneto- 
transport measurements have provided some support for 
the existence of two dimensional conduction channels 
originating from the surface states 0, 0]- Unfortunately, 
the transport measurements of TI have been masked by 
the conductivity of the bulk, leading to substantial diffi- 
culties in resolving the properties of the surface. 

In this paper we report on the response of Dirac 
fermions under crossed electric and magnetic fields. Un- 
like massive fermions, the conductivity, in a quantizing 
magnetic field, has an oscillatory dependence on the ap- 
plied electric field, wherein the even Landau levels have 
lower conductivity than their neighboring odd ones. This 
behavior is a consequence of the electric field dependence 
of the energies and wave-function of the Landau levels 
that arise due to the Lorentz invariance, allowing one 
to transform away the electric field in favor of an effec- 
tive magnetic field. The same property allows for novel 
features in the steady state thermopower unique to the 
linear dispersion. 

For linear spectrum tuning the in-plane electric field 
leads to a collapse of the Landau levels Q, (un) squeezing 
of the oscillator states and unusual dielectric breakdown. 
These phenomena are robust as long as the magnetic unit 
cell is much larger than the crystallographic unit cell but 
smaller than the system size. Here we focus on the behav- 
ior of the conductivity and thermopower of the 2D Dirac 
fermions, with spin orbit interactions, as function of in- 



plane electric field. For weak fields we use the standard 
linear response theory while steady state properties 
are computed using the density matrix approach [lOj . Our 
chief conclusion is the prediction of magneto-oscillations 
as a function of applied electric field at fixed chemical 
potential while no oscillations is expected at fixed parti- 
cle density. In the clean limit we find a correction to the 
universal values of thermopower, for chemical potential 
at the Landau level, whose field dependence provides a 
signature of the surface state. For weak spin orbit scat- 
tering, our results are also applicable to graphene. 

For an electric field applied along the surface of 3D 
TI and magnetic field perpendicular to the surface, the 
Hamiltonian is 



Ho = J d fxjr (?) vpa ■ U — jj,I — g^B® ■ B — e£x ip (r) 

where ip (r) is the two component electron wave func- 
tion, vp is the Fermi velocity, /i is the chemical potential, 
a = {o- x ,a y , a z } are the Pauli matrices representing spin, 
/ is a 2 x 2 unit matrix, II = —ihV — eA is the canonical 
momentum, g is the gyromagnetic ratio, /xs is the Bohr 
magneton, and £ is the magnitude of the electric field 
applied along the x axis. The magnetic field is pointing 
along z-direction, perpendicular to the surface, and we 
choose the gauge field A — (0, Bx, 0). The Lorentz in- 
variance of the Dirac equation implies that we can boost 
along y direction to a frame where the electric field does 
not appear in the Hamiltonian. The transformations is 
x' = x, y' — 7 (y + ftvpt) , t' = j(t + in spatial coor- 
dinate, £' = j(S - Pv F B),B' = 7(5 - p£/v F ) in fields, 
and xj} (r ) = cxp[a y tanh _1 (/3)/2]-0(r) in wavefunctions 
where 7 = l/^/l — (3 2 and (3 = £/vpB. In the boosted 
frame we solve for the eigenstates and transform back to 
get the exact eigenstates in finite electric field in the lab 
frame. These transformations are valid as long as fj ^ 1. 
For fj > 1 the transformations lead to scattering states. 
The transport properties for this case are beyond the 
scope of this investigation. In the lab frame the energy 
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eigenvalues are 



E n = sgn(n) . \/2 \n\ + n^B ~ f3hv F k y , n / 
7 ' "b 

= —■y^ 1 gp B B — fihvpky ,n = Q 



(1) 



with l B = \Jh/eB as the magnetic length, kq — 
g 2 p 2 B /hv 2 F e, and n as Landau level index for Dirac 
fermions. The corresponding spatial wave functions for 
ri 7^ in the lab frame are 



"0n (r) 



N n e 



zk y y 



cosh(|)<^|„| +za„sinh(|)(/)|, i |_ 1 
-isinh(|)^| n | +Q„cosh(|)0|„|_ 1 



(2) 



where we have used k y y — Et = k' y y' — E't' and take 
the spatial part of the wavef unction. N y is the nor- 
malization factor coming from periodic boundary con- 
dition in the y direction. We denote the parameters 

in the lab frame by N n = 1/yl + \a n \ 2 , tanh(6>) = 

/3, a n = —i (^sgn(n)^l + jc 2 B 2 /\n\ + ^/jcB / 



= Isg/J-B I '\2hvp , and 



-H\„\ (f n )e 



is given by 
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order Hermite polynomial, £ n = { 
with l n = sgn(n)pi B y/j^2 \n\ + n ~/B 
-~f3l B y/K J 2 B. For n = 



ui is the |n|th 
and Iq = 



ipa (r) = 



cosh(<9/2) 
-isinh(6»/2) 



ik y y 




(3) 



In the lab frame each component is a linear combination 
of Hermite polynomials of different index n. Moreover 
/3, being a function of the electric field, varies from zero 
to 1. Since 7 diverges as j3 approaches 1, the Landau 
levels will collapse Q. The fact that the relative strength 
of the orbital and Zeeman terms can be manipulated by 
the electric field or the magnetic field allows richer de- 
pendence of transport coefficients in the surface state of 
topological insulators. To compute conductivity we use 
the density matrix approach 10] to obtain the charge cur- 
rent. The electrical conductivity tensor a a p is given by 
J a = GapEp. For an ensemble described by the density 
matrix p, we get 



J(f) 



-evF 



mPn 



(4) 



To obtain transport coefficient we introduce a scattering 
potential V s (r). The density matrix satisfies 



H = H Q + V S 



(5) 



A Laplace transform yields —isP(s) — [P(s),H] — 
ip(0), where P(s) = / °° e~ st p(t)dt is the Laplace trans- 
form of p. We decompose the P{s) into three parts 



Pmn(s) = \(f n 5rnn + D mn (s) + G mn (s)) where f n /s 

is the Laplace transform of the initial density ma- 
trix which is assumed to be the Fermi Dirac distribu- 
tion /„ = f(u) n ) = l/(exp[(o; n - p)/k B T] + 1) with 
uj n = 8gn(n) ^)' 2 F lB a/2 \n\ + k Bj. D mn is the part of 

k„ of the 



the density matrix which is nonzero only if ny 
states m and n are identical and G m n is nonzero only 
if k y of the states m and n are different. Expand- 
ing to linear order in the scattering matrix the den- 
sity matrix is[HJ G mn = ^^-fThZk™-----™ and Dmn ~ 



v mv v v 



V, 



Ju 



(3hv pkyjnv —is ixi un — phvpk ylln —is ' H*^^ 

is the matrix element of the scattering potential, 

^mn ^nj ky nu ky n ky V and fav fa Jv 

For an arbitrary potential we write its Fourier com- 
ponent V s (r) = ~YliqVq e%qr an d the matrix element 
V nm = J dfipl(r)V (r)i> m (r) = J2q V qnm- In this paper 
we study the case of random impurity potential where 
V s (r) = J^i Vod(f— r\). For a random array of scatter- 
ers, the phase factors appearing in G nm will average to 
zero. We only need to consider matrix elements in the 
current operator where the two states have the same k y . 

We compute the steady state current by taking the 
long time limit or steady state density matrix p(i — > 
00). This is obtained by taking Laplace transform pa- 
rameter s through positive zero as p m n{t — >■ 00) = 
lim s _ ! .Q+ sP mn {s). For off diagonal elements of the den- 
sity matrix only D nm is relevant. This dissipation term 
is controlled by the imaginary part and the expression is 



D„ 



f3hVFk mpy ) 



(6) 



where k nmy is the difference in value of k y of the states n 
and m. Take /3 <C 1 in Eq.Q and Eq.© we obtain the 
diagonal conductivity a xx = J x /E in the linear response 
regime as 



^27r|K,| 2 (fcWB) ; 



,df 
duj 



S(u n u) (7) 



which is the same as the conductivity derived from lowest 
order Kubo formula with rij denoting impurity density. 
In Eq.([7]) we have used lowering operators on V n +i lV to 
connect it with V nu as done in Ref.flo| for 2DEGs case. 
For large in-plane field we compute numerically the di- 
agonal conductivity as a function of chemical potential 
in different in-plane electric field for both conventional 
2DEGs and Dirac fermions. 

The results for longitudinal conductivity for 2DEGs 
and Dirac fermions are shown in the left and right panel 
of figH] respectively. Within our approximation the Lan- 
dau levels have zero width and the conductivities show 
steps as a function of chemical potential (see the top left 
and right, and the middle left panels of figHJ. Note that 
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all features that we discuss have sharp transitions across 
Landau levels. This is related to the lack of impurity and 
thermal broadening. The results are valid as long as the 
broadening is less than the Landau level spacing. For 
2DEGs the conductivity either increases or decreases as 
a function of carrier density depending on whether e ^ B 
is 1 or < 1 . In the former case the increase of available 
phase space wins out against the decrease in the matrix 
elements V np as a function of chemical potential, while 
for the latter the opposite is true. The non-monotonic 
behavior is also expected as a function of electric field 
strength as seen in the bottom left panel. 

The behavior of the Dirac spectrum is far richer as 
there are a number of competing effects. The Fermi func- 
tion restricts contribution to the density matrices from 
levels near the chemical potential. Since the levels get 
closer for large n, the energy denominator in Eq.(|6]) de- 
creases, leading to an increase in the density matrix. The 
phase space satisfying the delta functions also increases 
due to the decrease in the level spacing at large n fur- 
ther enhancing D nm , On the other hand the matrix el- 
ements, V np , decrease for large chemical potentials. To 
understand why, note that the two states n and p differ 
in their k y values. Consequently their centers are sepa- 
rated in real space. As the chemical potential increases, 
the difference in k y values and hence the distance in real 
space, goes down. The two states are orthogonal, if they 
had the same center and thus the decrease at large chem- 
ical potentials. There are two further competing effects 
in the presence of an electric field. An increase in the dis- 
tance between the centers due to the change in l n com- 
petes with the decrease in the difference in k y due to the 
Landau level collapse. The net result of all these vari- 
ations is an even-odd effect, where the conductivity of 
the odd Landau levels are larger than their neighboring 
even levels as seen in the top right panel of fig[T]ll]. At 
large chemical potentials, the effect is washed out as the 
scattering process connects multiple levels with different 
symmetries. 

As evident in fig [1] and Eq.©, the order of magnitude 
of the conductivity is set by the impurity potential and 
impurity density within the perturbative approach em- 
ployed here. To get to conductivities of e 2 /h, a larger 
value of the impurity potential is required, beyond the 
scope of this calculation. Nevertheless, the Landau level 
collapse and the even-odd effect are more general prop- 
erties arising from the linear dispersion. To what extent 
these results will be modified within a fully self-consistent 
treatment is a subject of future investigation. 

For Dirac fermions the conductivity also shows oscil- 
latory behavior for fixed chemical potential as shown in 
middle right of figfT] This quantum oscillation has to do 
with the Landau level collapsing |8| with increasing field 
mainly due to the 7 -3 / 2 factor in u) n . The even-odd ef- 
fect translates to an oscillation as opposed to a series of 
plateaus. For fixed charge density, as implemented in the 
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FIG. 1: (Color online) Top: a xx v.s. fi/hu c . Top right: the 
Dirac fermions with E/vfB — 0.4 (blue circle) and 0.6 (pur- 
ple square). Top left: 2DEGs with eEl B /hw c ~ 0.77 (blue 
circle), 0.88 (purple square), and 0.99 (brown diamond). Mid- 
dle right: a xx v.s. E for chemical potential fi ~ 1.7/ftu c (blue 
circle), 2.2 (purple square), and 2.8 (brown diamond). The 
Landau level positions are shown explicitly for /i ~ 1.7/Hu) c - 
Middle left: a xx v.s. /i/ftu c for 2DEGs with eEl B /huj c ~ 44 
(blue circle), 66 (purple square), and 77 (brown diamond). 
Bottom right: a xx v.s. E(W 5 V/m) for fx = W2 (blue cir- 
cle), n = W4 (purple square), and fi = Lie (brown diamond). 
Bottom left: a xx v.s. E(W 5 V/m) for 2DEGs with chem- 
ical potential fj, = 4.5ftu c (blue circle), /j, = 5.5hu c (pur- 
ple square), and /i = 6.5ftu c (brown diamond). The order 
of magnitude ( 10" 18 and 10" 21 for Dirac and 2DEGs re- 
spectively) is set by the scattering matrix element given by 
Env \Vn*\ 2 6(E n - fi)/huj c with V ~ 0.18fc c at T = IK for 



the Dirac spectrum and Vo 
2DEGs case. 



Q.llfiwc at T = IQ-'K for the 



experiment by fixing the gate voltage, no oscillations are 
expected (bottom right panel of figO} as the degeneracy 
of the Landau level is independent of the applied electric 
field. 

Another interesting measurable quantity is the steady 
state thcrmopower. Consider a system with finite current 
due to an applied voltage V. We impose small source- 
drain temperature difference, such that the temperature 
gradient can be treated in perturbation, with a mean 
temperature T. The particle current density Jf is(l2| 



P(l) 



J, 



P(2) 



with J, 

7^(2) 



p(i) 



— J l /e where J 1 is the electric current, and J i is the 
particle current coming from temperature gradient. The 
electric conductivity a-ij = —{e 2 /T)L\j. Note that Ljj 
is not given by the Kubo formula as — W is not small 
enough to be in the linear response regime. We define 
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the out of equilibrium steady state thermopower Sij as 
the response to a small temperature gradient and small 
perturbation VV around the applied VV. The out of 
equilibrium thermopower is 



(8) 



Since the temperature gradient is treated within linear re- 
sponse and the large uniform electric field does not break 
time reversal symmetry, the L\j is L\ 2 — L\j + with 
L\] = linw,^ / °° dte^Tr(p[Jf(t),j9( )]) repre- 
senting the Kubo term and Ptj = ^Tr(p{fj Jf\) the 
correction term due to diamagnetic currents [12j. Here 
jf = a % ^) is the particle current operator, fi is the 
position operator, and jf = (Hjf + JfH)/2 — pjf is 
the heat current operator, p is the steady state density 
matrix. The time dependence is expressed in the Heisen- 
berg representation. The position operator is related to 
the particle current operator as Jf = -^-[H,fi]. In the 
clean limit where impurity potential energy V — > 0, the 
diagonal conductivity satisfy a xx = a yy = even in the 
presence of finite electric field (with E < vpB). Within 
this clean limit the low temperature thermopower in the 
finite electric field for 2D Dirac spectrum is shown in 
figHJ With increasing electric field we see that the ther- 
mopower increases and begins to deviates from its uni- 
versal value when f3 ~ 1CP 3 (E ~ 10 3 V/m), which is 
much smaller than the linear response regime seen in the 
transport measurement of current as shown in the bot- 
tom right of figHJ 

In conventional 2DEGs the thermopower (following 
the derivation in [10|]) at finite electric field, low tem- 
perature and for p ~ fiw c {N + 1/2) = En is S xx ~ 

-^rfy - ^iMtt) • The second tcrm kads t0 

an increase in the large electric field. In the linear re- 
sponse regime the 2DEGs has universal thermopower 
S xx ~ — TTfxr; at p ~ En, independent of the B field 
strength. From Eq.© the thermopower in the linear 
response regime for clean 2D Dirac gas with negligible 
-2 



Zeeman term is S r 



eT 



tanhl 



which is the 



same as derived by Kubo formula in At p ~ ui n 
the thermopower S TT ~ - 8 is universal in the lin- 

1 ^ en 

ear response regime. For chemical potential in between 
Landau levels the lack of density of states leads to vanish- 
ing thermopower. This is related to the lack of impurity 
and thermal broadening of the Landau levels. For suffi- 
cient large in-plane field the universal peaks in the ther- 
mopower are overwhelmed by the entropic contribution 
from the field, similar to its counterpart in the 2DEGs 
where the E 2 correction gives rise to overall increase in 
the thermopower. For materials with large gyromagnetic 
ratio or large Zeeman effect the peak positions in S xx 
depends on the strength of magnetic field and the peak 
height is not universal Q. As the nonequilibrium feature 




FIG. 2: (Color online) Negative thermopower —S X x(kB /e) 
v.s. chemical potential p/ftw c for the clean 2D Dirac fermions. 
Temperature is fixed at 10 _2 &j c /fcs. g = 2 and E = 10 2 V/m 
for blue circle line, E = l0 3 V/m for purple square line, E = 
10 4 V/m for brown diamond line. 



depends on the ratio between electric and magnetic field, 
larger magnetic field results show features similar to the 
equilibrium system. 

In summary, we use Lorentz boost to find the exact 
eigenstates in the presence of crossed electric and mag- 
netic field for Dirac fermions realized on the clean sur- 
face of 3D TI. The effect of weak random impurities are 
considered perturbatively through density matrix formu- 
lation. In the linear response regime the conductivity ob- 
tained is the same as the lowest order Kubo formula. For 
general in-plane field strength we carry out the numerical 
computation for both Dirac fermions and 2DEGs. Quan- 
tum oscillation is observed as a function of in-plane field 
for fixed chemical potential due to the Landau level col- 
lapse, which serves as unique signature of Dirac fermions. 
In addition to the charge transport, we define the steady 
state thermopower to extend the thermoelectric measure- 
ment to the nonequilibrium steady state. The external 
field, acting as the driving source for entropy production, 
enhances the overall thermopower. For sufficient large 
field the universal peaks in the thermopower are wiped 
out. 

Finally we discuss the observability of the effects dis- 
cussed. While the even odd effect is visible for clean 
samples (i.e. level broadening much less than level spac- 
ing and well defined parities of the wave- functions), a full 
self consistent calculation is currently being pursued to 
determine the regime of breakdown 15 1 and will be 
published elsewhere. Such a study will also help estab- 
lish the extent of the validity of Lorentz transformation 
for systems with dissipation. Joule heating can poten- 
tially wipe out the effect. The net heat flux produced 



by a lOyum x lOyum device is O.lmW at E ~ 10 5 V/m 
(a ~ e 2 /h). As long as the cooling is efficient enough 
to overcome this, the electric field effects are observ- 
able. Magneto-oscillations in linear response have al- 
ready been measured both in topological insulators @ 
and graphenejlij]. The latter has low enough minimum 
conductivity to make these measurements feasible. 
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